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ABSTRACT 

SMATASY is an interface for the zpl Tj^ji package and may be used for the model independent 
description of the Z resonance at LEP 1 and SLC It allows the determination of the Z mass 
and width and its resonance shape parameters r and j for cross-sections and their asymmetries. 
The r describes the peak height and j the interference of the Z resonance with photon exchange 
in each scattering channel and for ax-, crpB, crir, <^poi etc. separately. Alternatively, the helicity 
amplitudes for a given scattering channel may be determined. We compare our formalism with 
other model independent approaches. The model independent treatment of QED corrections 
in SMATASY is applicable also far away from the Z peak. 



^ Present address: Particle Physics Division, CERN. 



PROGRAM SUMMARY 



Title of program: SMATASY version 2.2 
Catalogue number: 

Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland (see 
application form in this issue) 

Licensing provisions: none 

Computer for which the program is designed: any computer with a FORTRAN 77 compiler 

Operating system under which the program has been tested: UNIX 

Programming language used: FORTRAN 77 

Memory required to execute with typical data: 34000 words 

No. of bits in a word: 32 

No. of lines in distributed program, including test data, etc: 1204 

Keywords: electron-positron annihilation, scattering matrix, cross sections, asymmetries, model 
independent ansatz, radiative corrections 

Nature of the physical problem: The program determines cross sections and asymmetries in the 
vicinity of the Z resonance taking into account radiative corrections up to O(a^) . This 
model independent determination is based on the S-matrix approach. 

Methode of solution: SMATASY is designed as a new interface of the ZFITTER package . 
Therefore it relies on the same treatment of radiative corrections as ZFITTER . 

LONG WRITEUP 
1 Introduction 

With the rising precision of the experimental study of the Z resonance at LEP 1 and SLC, 
an interpretation of the data in a rigorously model independent way becomes more and more 
important. A wrong theoretical description of the data may lead to systematic non-observed 
shifts of the measured parameters which describe the Z. 

SMATASY is a program for the model independent study of fermion pair production at the 
Z resonance: 

e+e-^(7,Z)^/+/-(n7). (1) 

Under the assumption that QED and QCD are well understood theories, the total cross-section 
around the Z peak may be characterized by four real parameters: the mass Mz, width Tz, the 
residuum r of the Z resonance, and the strength of the 7Z-interference j. The parameters r and 
j are related to helicity amplitudes Rz'^ for e~^e~ annihilation into fermion pairs via Z-exchange. 
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Further, the measured cross-section includes photonic virtual and real bremsstrahlung correc- 
tions which may be described by a flux function Pt{s' / s) |l|]: 

4 ^ rds' 
<^t[s) = -na / — 

6 J S 

As mentioned, the photon exchange parameter r'' is assumed to be known, and s depends on 
the beam energy, s = AE"^. Further, at LEP 1 and SLC very precise measurements of various 
cross-section asymmetries are performed. An example is shown in figure |l|. In the vicinity of 
the Z peak, these asymmetries behave relatively smoothly and may be described by a simple, 
universal formula |@]: 

A{s) = Ao + C{s)A,(^-^-iy (3) 

The QED corrections are contained in the factor C{s). 

The model independent description of cross-sections around the Z resonance with account 
of QED corrections may be done with the zfI Tj^j^ branch ZUSMAT. Here, we describe 

the Fortran program SMATASY which is designed as an interface to zfI Tj^ji for the model 

independent description of asymmetries. SMATASY provides the full functionality of zpl Tj^ji 
with all its possibilities of flag settings, different treatments of photonic bremsstrahlung and 
QCD corrections. 

SMATASY is devoted to the following tasks: 

• determination of the Z-exchange parameters r and 7Z-interference parameters j for cross- 
sections and asymmetries, 

• determination of the asymmetry parameters Aq and Ai, 

• determination of the couplings Rz'^ of helicity amplitudes describing the Z-exchange ma- 
trix element, 

• model independent treatment of QED corrections for cross-sections and asymmetries. 

The basic formulae are introduced in section 2 and related to other approaches in section 3. 
The structure of the SMATASY package is explained in section 4 while the procedures are 
described in section 5. Appendix A contains a sample output of SMATASY. 
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sr + {s- Ml)j 



r 

~^ {s-MlY + MlTl^ 




(2) 



2 Basic formulae 

2.1 Cross-sections at the Z resonance 

The matrix element for the production of a fermion pair near the Z resonance may be described 
by four helicity amplitudes. In full generality, they have the following form |^: 

Mns) = ^ + ^ + t^i-'-r] ^ . = 1,...,4. (4) 
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Figure 1: The forward-backward asymmetry for the process e~^e fi~^fi near the Z peak. 

The position of the Z pole in the complex s plane is given by sz'- 

sz = mI- iMzTz. (5) 

The Rlj, and Rz'^ are complex residua of the photon and the Z boson, respectively. In Born 
approximation, they are real numbers. In the Standard Model, contributions from higher order 
corrections are incorporated and the residua become complex numbers. The coefficients 
in describe nonresonant contributions to the scattering process. In the Standard Model they 
arise from higher order corrections, e.g. from ZZ and WW box diagrams. The s dependence of 
the virtual electroweak corrections also contributes to them. 

To a very good approximation, the Taylor series in @j may be neglected. Its first term Fq* 
is constant. In a Z line shape fit, it will be strongly correlated with the photon exchange term 

r;/s ^ r;/m|, 

R^ = FA(s)|QeQf|. (6) 

The Qf is the electric charge of the final state fermions and Fa{s) the vacuum polarization of 
the photon: 

f,(M|) = ^^ = HL£5_, 0.0188. (7) 

a 128.86 ^ ' 

In practice it seems to be impossible to disentangle Fa and the nonresonating quantum cor- 
rection Fg . The next-to-leading background term is proportional to cr = (s — Sz)/M2 
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2{^/Mz - 1). At LEP 1 and SLC the bulk of data is taken at \^ - M z\ < Tz and a 
becomes less than 5 %. Having in mind that the Fj^^ are quantum corrections and proportional 
to a/vr, it is for all practical needs: 



s 
s 



s-sz Mz 



s - Sz' 



There are four residua Rz^ for the four independent helicity amplitudes in the case of massless 
external fermions: 



(9) 



The amplitudes A4^\s) give rise to four cross-sections which add up incoherently to the 
following measurable cross-sections: 
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(10) 



0"2 + CTs. 

Here, the 0"° is the total cross-section, a^^ defines the forward-backward asymmetry, cTp^j the 
final state polarization, ap^ the forward-backward asymmetry of the final state polarization 
etcR 



All these cross-sections may be parameterized by the following master formula pi] 



-IT a 
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The r]/ is the photon exchange term. 



A= T, FB, pol, FB-pol, . . . 



(11) 



.7/ 



rC/E{±l} 



i=0 



R 



QCD-, 



(12) 



vanishes for all asymmetric cross-sections, cy = 1,3 for leptons and quarks, respectively. QCD 
corrections for quarks are taken into account by the factor Rqcd of Hi- The Z-exchange 

^The agreements between several pairs of asymmetries will be disturbed after inclusion of QED corrections, 
see below. 
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residuum and the 7Z-interference are: 



3 a 



I ^ i=0 



(13) 



The factors {±1} in ([T^) and (|1^) indicate that the signs of Rl^ 
to the signs of ai in (plQl). 



R7 



and of Rz^ correspond 



2.2 QED corrections for cross-sections 



rp 

For the calculation of QED corrections in SMATASY the zpl Tj^ji environment is used. This 



is done by convoluting (|TT]) with radiator functions for initial and final state radiation and 
their interference. The initial and final state corrections with soft photon exponentiation to the 
cross-sections ctt, (Tpoh ^ir, (^ir-poi are described by: 



f ds' f s'\ 

J — a^{s')pT ( — j A = T, pol, Ir, Ir-pol. 



(14) 



Analogously, for the forward-backward difference cross-sections, 

a = FB, FB-pol, FB-lr. 



(15) 



We stress that Pfb{s' / s) ^ pt{s'/s). The contributions from initial final state interference 
bremsstrahlung are slightly more complex At the Z resonance they are numerically tiny 

as long as no strong cuts are applied. They will be neglected in the following. With QED 
corrections, the master formula may be rewritten as follows H: 
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The barred parameters contain correction factors with QED corrections: 
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where 



I 
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is - mIy + mItI 



f2 t^2 



18) 



m|)o 



X[l]. 



Here, the definition 



(19) 



is used. The QED correction factors are completely model independent, i.e. independent of the 
underlying dynamics of the scattering process. They depend on mass and width of the Z and 
on the handling of the photonic phase space, the inclusion of higher orders, and on acceptance 
cuts. The reader may wonder that the corrections C\ seem to be singular at ^/s = M z- This 
is not the case for the products C^(s)(s — Mz) which are physically relevant. As may be seen 
from the corresponding definitions, these products remain small (but potentially non- vanishing) 
when a/s approaches Mz- There the QED corrected cross-sections may be defined as (smooth) 
limits from the neighboring energies. 

The correction factors are shown in figures ^|-|^ without and with two different cuts on 
the photon phase space. All QED corrections are smooth and, with one exception, rather 
independent of s. Those to the Z exchange, C^, develop the radiative tail at the right hand 
side of the peak which at some value of s gets suppressed by the cuts. It may be further seen that 
the corrections to total cross-sections and those to forward-backward differences are not equal, 
although of similar size near the peak. They deviate more when more hard bremsstrahlung is 
possible 1^ . This explains the rise of their difference with the tail and the subsequent vanishing 
of it after the cuts become influential. 



2.3 Asymmetries around the Z resonance 

Without QED corrections, asymmetries are defined by: 



A=^T. 



(20) 



These asymmetries take an extremely simple form around the Z resonance 



1 + A 
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A^ + At 



Mr, 



(21) 
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Figure 2: The model independent QED correction factors without cuts. 
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Figure 3: The model independent QED correction factors with a cut on the energy of the 
bremsstrahlung photon. 
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Figure 4: The model independent QED correction factors with cuts on the acoUinearity ^ of 
the final state fermions and on their energy Ef. 
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The higher order terms may be safely neglected since {s/M^ — 1)^ 
coefficients have a quite simple form: 
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(22) 
(23) 



Here, the r|4^ is neglected in both Aq and Ai. Further, the definition 7^ 
is used. 



T^2 



0.75x10 
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2.4 QED corrections for asymmetries 

A typical cross-section asymmetry with QED corrections is shown in figure The Born asym- 
metry is linear around the Z resonance. QED corrections lead to deviations from this; specially 
pronounced at the right hand side of the peak. Although asymmetries get much smaller QED 
corrections than the cross-sections themselves, an analysis of the data would not be consistent 
without their correct treatment. 

With QED corrections, the master formula for asymmetries becomes: 



Aa{s) = A^ + At^-l\. (24) 




The coefficients An may be obtained from the An defined in section p.3| by replacing in their 
definitions the unbarred variables by barred ones. 

For the peak contributions to the forward-backward asymmetries the explicit expressions 

are: 



^0 



f 



^FB 'a 



CJp 4 + 0.001 

/ / 

^ 0.998^^^^^^ ^ a = FB, FB-pol, FB-lr. (25) 

4 + 0.001 4 

Aiming at an experimental accuracy of AA% = 0.001 at LEP 1, the deviation of C^b/Ct 
from unity has to be taken into account. Compared to Afb, Aps-poi and Aps-ir the expres- 
sions for Apoi, Air and Air-poi are simpler since the QED corrections to the numerator and the 
denominator are both of the total cross-section type. The leading term is: 

A 



° 4 + [C?/C^]724^ 



4 + 0.001 



A = pol, Ir, Ir-pol. (26) 
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The coefficients are in reasonable approximation: 
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(27) 
(28) 



The coefficients with a pronounced s dependence are 



C{s) ^ Oris) = ^, 



c 



FBK.S) 



FB 
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(29) 



FB 



The behaviour of Ct and C^s as functions of s and the dependence on cuts is shown in 
figure^. As mentioned aheady, the Ct(s/M^ — 1) does not vanish at ^/s = Mz but is extremely 
small. Aiming at an accuracy of several per mill, one may neglect the difference between Ct 
and CpB in the vicinity of the peak. The initial state QED corrections to the Z-exchange 
cross-section develop a radiative tail while those to the 7Z-interference do not. Due to this, 
their ratio C{s) gets damped at the right hand side of the Z peak. This damping was seen in 
figure 1. It may be assigned to QED corrections completely. At ^/s > Mz, the radiative tail 
may be avoided by a cut on the energy of the emitted photons: 



< A = 1 - 



Ml 



beam 



(30) 



At LEP 1 and SLC, where ^/s ~ Mz, this condition is rather restrictive; e.g. at ^/s = Mz+^Tz, 
it is A = 0.1. In figures 1 and 5, photons are cut with energies larger than 6 GeV. The radiative 

tail is suppressed by this if y/s' > Mz is ensured. At ^/s > Mz/ \/l ~ -E-^/i^beam = 97.5 
GeV this is the case. In the immediate vicinity of the peak one unavoidably measures data 
which contain radiative corrections. As may be seen from the figure, the other cuts (on the 
acollinearity and energy of the final state fermions) are similar P]. 



3 Other model independent approaches 
3.1 The BCMS approach 

In ^] the following model independent formula for the total cross-section has been advocated: 



\s-s^\- [My M-'z ' 'Mz 

The free parameters of the above expression may be calculated within the Standard Model with 
account of electroweak radiative corrections 
dence to the parameters in (|Tl]) [||: 



If— + . 
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QED- 



(31) 



T0|| . There is a simple one-to-one correspon- 



a2 
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Figure 5: The model independent QED correction to the slope of asymmetries near the Z 
peak. 
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Further, 



So = M| - ^ 



(33) 



The relation between the definitions of Sq and sz is given to a very good approximation by the 



following equalities, which also affect the coupling strength JIT 



Mz 



[1 + {T^/Mzf]-"^Mz 
[l + (rz/Mz)2]-5rz 
G^/(l + irz/MZ). 



Mz - 34 MeV, 
Tz - 1 MeV, 



(34) 



3.2 The OPAL approach 



A completely different approach has been chosen in |T2|. With an ad hoc ansatz, the effective 
couplings of the differential cross-section have been allowed to deviate from what is expected 
in the Standard Model in the approximation of effective fermion couplings: 
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Here, the definitions 
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(39) 
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are used for the effective couplings in the Standard Model. The parameters n in the definitions of 
the parameters C allow for deviations of the cross-sections and asymmetries from the Standard 
Model predictions where they are equal to one. 

The OPAL approach uses the zfI Tj^ji environment with minor modifications. 

Up to tiny terms which are presumably much smaller than the experimental accuracy (e.g. 
the parameter in the S-matrix approach), the OPAL approach is equivalent to the one 
advocated here. For the total cross-section and the forward-backward asymmetry the relations 
are: 
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(40) 



^cFa + ^'^mFA 



The additional factors of | in the last two relations are due to the different angular integrations 
over even or odd integrands with respect to cosd. 



3.3 Effective couplings 

Now we define the relation of the model independent approach of SMATASY to the use of 
effective weak neutral fermion couplings. The latter is realized in the zfI Tj^^ branch ZUXSA. 

In a simple quantum mechanical interpretation or an approximate quantum field theoretical 
one, the (complex) residua of the helicity amplitudes Rz^ may be expressed in terms of effective 
vector and axial vector couplings of the Z boson to fermions (which are basically real numbers) : 
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with the K of ( P7| ) and the couplings gj^ and cjy of (^). The parameters r^^ and Ja of (0) are: 
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(42) 
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and 
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(43) 



jpB-pol - ~ '^l^9\g{f{C Re + C im) ■ 



Additional factors of (3/4)''=^ are again due to the different angular integrations for contributions 
which are even or odd in cos "d. 

The 7Z- interference is proportional to C/je, while Cim are small corrections to it and to the 
resonance peak parameter: 



Cim = I QeQf I — 5>mFA(s), 
IVlz 

CRe = \QeQf\ ^eFAis). 



(44) 



3.4 Model independent QED corrections 

The QED correction factors in SMATASY are universal in the sense that they may be used also 
at energies far away from the Z peak and in other approaches. We give here one instructive 



example for this. In [jl3l, the QED corrections to the total cross-section o"t(s) and the integrated 
forward-backward asymmetry have been calculated analytically without cutsQ to order 0{a) 
for reaction (p. The Z exchange cross-section contribution has been presented there as follows: 



0-7 



\x{s)f{vl + al) [vj + a 



i + -(Qlno + QeQf ^4 + ^2 



(45) 



In the simplest case (no cuts, neglect of final state masses, no higher order corrections), the 
7^2 e.g. is the well-known QED final state correction |. The TYq contains the initial state 
corrections and 7^4 those from the initial final state interference. Similar but considerably more 
involved analytic expressions were derived for the forward-backward asymmetry. The following 
relation holds: 



CUs) 



a 

i + -n 

TT 



s) + ..., 



(46) 



where the dots stand for higher order corrections and a potential inclusion of final state radiation 
in the C^ and the for potential cuts. Following this example, the interested reader may 
find analogue relations for the other QED corrections. 



■^Thc corresponding analytic formulae with a cut on the energy of the bremsstrahlung photon may be found 
in the unpublished Fortran program ZBIZON. 
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4 Structure of the package 



rp 

For the installation of SMATASY the user has to replace subroutine BORN of zpl Tj^ji with 
subroutine BORN of SMATASY. 

To run SMATASY one has to initialize first zpl Tj^^ following the procedure described 
in P], section 6. Then, SMATASY is initialized by a call to subroutine ASYINIT. Subroutine 
ASYTEST illustrates the initialization procedure and performs a comparison of SMATASY with 
the other model independent approaches of zpl Tj^j^. The results are listed in an sample output 



in appendix |5.2.5| . 



The SMATASY package contains the following interface routines: 

• SMATASY - calculates total cross-sections and asymmetries as functions of the center-of- 
mass energy, the Z mass and width, the Z- and 7-exchange terms and the 7Z-interference 
terms; 

• SMATRZ - calculates total cross-sections and asymmetries as functions of the center-of-mass 
energy, the Z mass and width and the helicity amplitudes; 

• SMATAOl - calculates asymmetries as functions of the center-of-mass energy, the Z mass 
and width, the Z- and 7-exchange terms and the 7Z-interference terms for the total 
cross-sections and the asymmetry parameters, Aq and A^. 

Utility routines of interest for the user are: 

• CORQED - calculates the QED correction factors as functions of the center-of-mass energy 
and the Z mass and width; 

• RZFRVA - calculates the residua of the helicity amplitudes as functions of the Z mass and 
the effective couplings; 

• RJFRRZ - calculates the Z- and 7-exchange terms and the 7Z-interference terms as func- 
tions of the Z mass and width, the helicity amplitudes and the vacuum polarization; 

• AOIFRRJ - calculates the asymmetry parameters, Aq and A^, as functions of the Z- and 
7-exchange terms and the 7Z-interference terms; 

• ASYTRAF - performs the transformation of the Z mass, width and Fermi's constant between 
the two definitions in (041). 



Although our model independent ansatz implicitly assumes massless fermions since it is 
based on four different helicity amplitudes, corrections due to final fermion masses are applied 

rp 

in the sample output in order to be compatible with zpl Tj^^. However, the corrections for 

rp 

leptons and light quarks may be switched off by the zpl Tj^^ flag POWR. The deviations between 

rp 

different branches of zpl Tj^^^ itself and of the interface SMATASY are at most of the order of a 
few tenth of a percent. The most accurate asymmetry measurement at LEP 1 is expected for 
the forward backward asymmetry for leptons at the peak where a theoretical accuracy of less 
than 0.1 % is demanded. The internal deviation between different descriptions in the sample 
output for this quantity is about 0.01 %. 
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final state 


vv e"'"e T^T 


uu dd cc 


ss tt bb 




Bhabha 


INDF 


12 3 


4 5 6 


7 8 9 


10 


11 



Table 1: zfI Tj^ji convention for final state labels. 

5 Description of the procedures 

If not stated differently, the input and output arguments of the following subroutines are of the 
DOUBLE PRECISION type. 

5.1 Interface Routines of SMATASY 
5.1.1 Subroutine SMATASY 

Subroutine SMATASY is used to calculate total cross-section and asymmetries as functions of ^/s, 
Mz, Tz, r^, and r]/ . The first three coefficients of the Taylor series r^'"^ are also consid- 
ered. This refers to the cross-section parameterization ([Tl|) . For total cross-sections SMATASY 
is equivalent to the interface ZUSMAT, but note the different definition of j^. 

CALL SMATASY ( INDF , SS , SZMASS , SGAMZ , RT , JT , GT , FT , RA , JA , FA , I A , ASY* ) 

Input Parameters: 

INDF is the fermion index (see Table HD (INTEGER). 

SS is the center-of-mass energy, ^/s, in GeV. 

SZMASS is the Z mass, M^, in GeV. 

SGAMZ is the total Z width, Tz, in GeV. 

RT is the Z-exchange term for the total cross-section, r:^. 

JT is the 7Z-interference term for the total cross- sect ion, jj.. 

GT is the 7-exchange term for the total cross-section, rj/. 

FT is a vector of the first three coefficients , describing nonresonant contributions. 
RA, JA, FA are corresponding parameters for a cross-section difference a\. 
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lA defines which cross-section or asymmetry is calculated (INTEGER) 0: 



lA 



ITOT 


= - 


ax 


IFB 


= 1 - 


-Afb 


IPOL 


= 2 - 


~^ •^pol 


IFBPOL 


= 3 - 


ApB-pc 


ILR 


= 4 - 




IFBLR 


= 5 - 


ApB-lr 


ILRPOL 


= 6 - 


~^ A-lr-pol 



Output Parameter: 

ASY is the total cross-section or an asymmetry selected by lA. 

5.1.2 Subroutine SMATRZ 

Subroutine SMATRZ is used to calculate total cross-sections and asymmetries as functions of ^/s, 
M z, T z and the residua of the helicity amplitudes, Rz^, as introduced in (||). The nonresonant 
contributions, F|j in (|) are set equal to zero. According to (H), instead of the vacuum 
polarization, F4, is used as free parameter. 



CALL SMATRZ (INDF,SS,SZMASS,SGAMZ,RZ0,RZ1 ,RZ2,RZ3,FA,IA,ASY*) 



Input Parameters: 

INDF, SS, SZMASS, SGAMZ, lA have the same meaning as in subroutine SMATASY, ex- 



plained in section 5.1.1 



RZO, RZl, RZ2, RZ3 are the residua of the hehcity amplitudes, Rz^ (C0MPLEX*16). 
FA is the vacuum polarization, Fa{s) (C0MPLEX*16). 
Output Parameter: 

ASY has the same meaning as in subroutine SMATASY. 

5.1.3 Subroutine SMATAOl 

Subroutine SMATAOl is used to calculate asymmetries as functions of y^, M^, Tz r^, jx, Tj/ 
and the asymmetry parameters, Aq and A^, introduced in |23|). The nonresonant contri- 
butions are neglected. 



CALL SMATAOl (INDF , SS , SZMASS , SGAMZ ,RT , JT, GT, AO , Al , lA , ASY*) 



Input Parameters: 

INDF, SS, SZMASS, SGAMZ, RT, JT, GT, lA have the same meaning as in subroutine 
SMATASY, explained in section |5.1.1| . 



''The reserved FORTRAN constants (e.g. ITOT) for the different possible lA values are given in the second 
column below. 

'^An asterisk (*) following an argument in a calling sequence is used to denote an output argument. 
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AO, Al are the asymmetry parameters Aq and A"^. 
Output Parameter: 

ASY has the same meaning as in subroutine SMATASY. 

5.2 Utility Routines of SMATASY 

5.2.1 Subroutine CORQED 

Subroutine CORQED calculates Cj., C^, C\ and as functions of a/s, Mz and Tz, according 
to flTT' 



CALL CORQED (INDF, SS , SZMASS , SGAMZ , CAR* , CAJ* , CAG* , CAO* , lA) 



Input Parameters: 

INDF, SS, SZMASS, SGAMZ and lA have same meaning as the parameters explained in 



Section E7L 



Output Parameters: 

CAR is the QED correction factor, C^, for the Z-exchange term, r{^, selected by lA. 

CAJ is the QED correction factor, C\, for the 7Z-interference term, j^, selected by lA. 

CAG is the QED correction factor, C^, for the 7-exchange term, r]/ , selected by I A ^ 

CAO is the QED correction factor, C^, for the nonresonant contribution, r^, selected by 
lA. 

5.2.2 Subroutine RZFRVA 

Subroutine RZFRVA is used to calculate Rz'^ as functions of Mz, G^, ^y, cjy and us- 
ing ( pTI) . Mz is related to Mz by (@). The subroutine cannot be used for the inclusive hadron 
channel (INDF=10). 



CALL RZFRVA (INDF, ZMASS , GMU, GVE, GAE, GVF, GAF , RZO* , RZl* ,RZ2* ,RZ*) 



Input Parameters: 



INDF corresponds to the parameter with the same name in Section ^ 
ZMASS is the mass of the Z boson, Mz- 
GMU is Fermi's Constant, G^. 

GVE is the effective vector coupling of the electron, ^fy. 
GAE is the effective axial vector coupling of the electron, 



Cl = for A^T. 
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GVF is the effective vector coupling of tfie final-state fermion, g^. 
GAF is the effective axial vector coupling of the final-state fermion, gj^. 
Output Parameters: 

RZO, RZl, RZ2, RZ3 correspond to the parameters explained in section [5.1.2| . 

5.2.3 Subroutine RJFRRZ 

Subroutine RJFRRZ is used to calculate r-^, and r]/ as functions of Mz, ^z, R-z^ and -Fa(s), 
according to (0). The subroutine cannot be used for the inclusive hadron channel (INDF=10). 



CALL RJFRRZ (INDF , SZMASS , SGAMZ , RZO , RZl , RZ2 , RZ3 , FA , RR* , J J* , GG* , lA) 



Input Parameters: 

INDF, SZMASS, SGAMZ, RZO, RZl, RZ2, RZ3, FA and lA have same meaning as the parameters 



explained in Section p3 . 
Output Parameters: 

RR is the Z-exchange term, r^, for the cross section, o"^, selected by lA. 
JJ is the 7Z-interference term, j^, for the cross-section, cr^, selected by lA. 
GG is the 7-exchange term, r]/ , for the cross-section, a a, selected by I A. 

5.2.4 Subroutine AOIFRRJ 

Subroutine AOIFRRJ is used to calculate Aq and as functions of Mz, Pz, rj,, j^, rjf, r^and 



jj^, according to ( p2 



CALL AOIFRRJ (INDF , SZMASS , SGAMZ ,RT, JT , GT, RA , JA , AO* ,A1*) 



Input Parameters: 

INDF, SZMASS, SGAMZ, RT, JT, GT, RA, and JA have same meaning as the parameters ex- 



plained in Section |5]1 
Output Parameters: 



AO and Al have the same meaning as the parameters explained in Section 5.1 



5.2.5 Subroutine ASYTRAF 

Subroutine ASYTRAF is used to perform the transformation of Mz, Pz, and from the nota- 
tions where the width of the Z is ^/s dependent to the parameters in S- matrix notation, Mz, 
Tz, and G^, according to 



CALL ASYTRAF (ZMASS , GAMZ , GMU, SZMASS* , SGAMZ* , SGMU*) 



Input Parameters: 
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ZMASS has the same meaning as the parameter in Section [5.2.2 
GAMZ is the total width of the Z boson, Tz- 
GMU is the Fermi Constant, G^. 
Output Parameters: 



SZMASS and SGAMZ have the same meaning as the parameters explained in Section |0 . 



SGMU is Fermi's Constant, (COMPLEX* 16). 
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TEST RUN OUTPUT 



** This 



is ZFITTER version 4.53 ** 

92/10/14 ** 



-k ^ 
-k -k 


The authors of 


the ZFITTER package are: 


* 


-k * 


D . Bardin 


(Dubna) 




■k -k 


M. Bilenky 


(Dubna) 




■k -k 


A . Chizhov 


(Dubna) 


-k 


■k -k 


A . OTshevsky 


(Dubna) 


-k 


■k -k 


S . Riemann 


( Zeuthen) 




* * 


T . Riemann 


( Zeuthen) 


* 


* * 


M. Sachwitz 


( Zeuthen) 


* 


* * 


A . Sazonov 


(Dubna) 


* 




Yu . Sedykh 


(Dubna) 


* 


* * 


I . Sheer 


(UC San Diego) 


* 
* 



** Questions and comments to ZFITTER@CERNVM. CERN . CH ** 
****************************************************** 



ZFITTER cut values: 

INDF ICUT ACOL EMIN S_PR ANGO ANGl 

-1 .00 .0000 .0000 .00 180.00 
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.00 



.0000 



.0000 



.00 180.00 



***************** 
******************** 

** This 

•kii'k'k'k'kii'k'k-k-k'k'k'k'k'kii 

** The authors of 

■A- *■ 

S.Kirsch 
** T. Riemann 

•kii'k'k'k'k'k'k'k-k'k'k'k'k'k'k'k 

** Questions and 
***************** 



************************************* 
************************************* 

is SMATASY version 2,1 
94/05/27 



* * 

************************************* 

* * 

* * 

* * 

* * 

* * 

************************************* 
comments to KIRSCHSCERNVM. CERN . CH ** 
RIEMANN@CERNVM.CERN.CH ** 
************************************* 



the SMATASY package are: 

(DESY IfH-Zeuthen) 
(DESY IfH-Zeuthen) 



ZFITTER 


flag 


values : 












AFBC : 


1 


ALPH 


ALST: 


1 


AMT4 : 


3 


BORN: 





BOXD : 





CONV 


FINR: 


1 


F0T2 : 


3 


GAMS : 


1 


INCL: 


1 


INTF 


BARB: 


2 


PART: 





POWR: 


1 


PRNT: 


1 


QCD3 


1 QCDC: 


1 


VPOL: 


3 


WEAK: 


1 



ZMASS = 



91.19500 



HMASS = 300.00000 



ALPHST 
QCDCOR 



.12000 
1.03954 



TMASS = 140.00000 

ALFAS = .12000 

SIN2TW = .22815 

QCDCOB = 1.04211 



Mass Z I Gamma Z | Gmu | 
91.1951 2 .487 I . 11664E-04 I 
91 .161 I 2.486 I . 11655E-04 | 



I S-Matrix | 
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SQRT (S) 
< 



90.1611 

Xs- 



-Asymmetry- 



DF 


1 ZUXSEC 


ZUXSA 


ZUSMAT 


SMATASY 1 


ZUXSA 


ZUTAU 




SMATASY 






















tot 1 


t JD 


pol 1 fbpol 1 


r JD 


pol 


fbpol 1 


Ir 


rJoir 


1 




8631 




8634 




8603 


.8603 1 


-.0961 




-.0963 






-.1191 


-.0892 


2 




8530 




8533 




8503 1 


. 8503 1 


-.0973 




-.0974 






-.1205 1 


-.0903 


3 




8465 




8467 




8451 1 


.8451 1 


-.0978 


-.12321 -.09201 


-.0979 


-.1212 


-.0907 1 


-.1212 1 


-.0907 


4 


1 2 


9770 


2 


9777 


2 


9715 1 


2.97151 


.0000 




.0000 






-.09931 


-.4686 


5 


1 3 


8128 


3 


8135 


3 


8076 1 


3.80761 


.0653 




.0653 






-.10991 


-.6713 


6 


1 2 


9645 


2 


9651 


2 


9516 1 


2 . 9616 1 


.0000 




. 0000 






-.0995 1 


-.4696 


7 


1 3 


8120 


3 


8127 


3 


8068 1 


3.80681 


.0653 




.0653 






-.1100 1 


-.6714 


9 


1 3 


7430 


3 


7436 


3 


7573 1 


3.7573 1 


.0650 




.0646 






-.1099 


-.6662 


10 


1 17 


3093 




0000 


17 


3048117.30481 


.0000 




. 0428 






-.1063 1 


-.6010 



Irpol 



, 1312 



SQRT 
INDF 

1 

2 
3 
4 

5 
6 
7 
9 
10 
SQRT 

INDF 

1 
2 

3 
4 
5 
6 
7 
9 
10 



(S) = 91 
I ZUXSEC I 



I 



I 6 



4854 I 
4754 I 
4676 I 
1892 I 
6679 I 
5 . 1737 I 
6.6671 I 
5492 I 
2471 I 
= 93, 



I 30 
(S) 

I ZUXSEC I 
I I 
. 6510 I 
. 6415 I 
. 6358 I 
2 .2371 I 
2 . 8552 I 
2 .2263 I 
2 . 8544 I 
2.8040 I 



1611 
— Xs — 
ZUXSA 

1.4858 
1 . 4759 



4681 
1904 
5691 
1749 
5583 
5504 
. 0000 
1611 
— Xs— 



I 12 . 9770 I 





>< 




Asymmetry- 










> 


ZUSMAT 1 SMATASY | 


ZUXSA 1 


ZUTAU 1 




SMATASY 






1 


1 


tot 1 


fb 1 


pol 1 fbpol 1 


fb 


pol 


fbpol 1 


Ir 1 


fblr 1 Irpol 1 


1.4808 1 


1.48081 


-.0056 1 


1 1 - 


.0056 






-.1293 1 


-.09681 1 


1.47091 


1 . 4709 1 


-.0057 1 


1 1 - 


.0056 






-.1301 1 


-.09741 1 


1.4653 1 


1.46531 


-.0057 1 


-.13261 -.09901 - 


.0056 


-.1305 


-.0977 i 


-.13051 


-.09771 -.00781 


5.17981 


5.17981 


.0527 1 


1 1 


.0527 






-.1273 1 


-.47341 1 


6. 6588 


6. 6588 1 


. 0858 1 


1 1 


. 0858 






-.1295 


-.67231 1 


5.1589 


5 .1589 1 


. 0528 1 


1 1 


. 0528 






-.1275 


-.47391 1 


6 . 6579 


6 . 6579 1 


.0858 1 




. 0858 






-.1295 


-.67241 1 


6.5735 1 


6.5735 1 


.0855 1 


1 1 


. 0851 






-.12951 


-.66711 1 


30.2387130.23871 


.0000 1 


1 1 


. 0743 






-.1288 1 


-.60341 1 



>< Asymmetry 



ZUXSA 


ZUSMAT 1 SMATASY | 


ZUXSA 








1 


tot 1 


fb 




6512 




6488 1 


. 6488 1 


.1141 




6417 




6394 1 


. 6394 1 


.1158 




5350 




5345 


. 6345 1 


. 1157 


2 


2377 


2 


2329 


2 . 2329 1 


. 1217 


2 


8557 


2 


8512 


2 . 8512 1 


. 1126 


2 


2258 


2 


2239 


2 . 2239 1 


. 1221 


2 


8550 


2 


8505 


2 . 8505 1 


. 1126 


2 


8045 


2 


8138 1 


2 . 8138 1 


.1125 




0000 


12 


9720 1 12 . 9720 | 


.0000 



ZUTAU 



pol 



1416 



I 



fbpol I 
I 



I 



.1055 



fb 

. 1145 
. 1162 
. 1169 
. 1218 
. 1126 
. 1222 
. 1127 
. 1120 
. 1157 



SMATASY 
pol 



,1396 



fbpol 



-.1043 



Ir 1 


fblr 1 


-.1367 1 


-.1021 1 


-.1387 1 


-.1036 1 


-.1395 


-.1043 1 


-.1525 


-.4724 1 


-.1552 


-.6705 1 


-.1532 


-.4738 1 


-.1553 


-.6707 1 


-.1555 1 


-.6658 1 


-.15791 


-.6018 1 



Irpol 



> 

I 
I 
I 
I 

1557 I 
I 
I 
I 
I 
I 
I 



